(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 


(19) World Intellectual Property Organization 
International Bureau 



(43) International Publication Date (10) International Publication Number 

10 January 2002 (10.01.2002) PCT WO 02/03041 Al 


(51) International Patent Classification 7 : GO 1H 1/00 

(21) International Application Number: PCT/GBO 1/03020 

(22) International Filing Date: 5 July 2001 (05.07.2001) 

(25) Filing Language: English 

(26) Publication Language: English 


(30) Priority Data: 

0016561.3 


5 July 2000 (05.07.2000) GB 


(71) Applicant (for all designated States except US): ROLLS- 
ROYCE PLC [GB/GB]; 65 Buckingham Gate, London, 
Greater London SW1E 6AT (GB). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): ANUZIS, Paul 

[GB/GB]; 86 Locko Road,: Spondon, Derby, Derbyshire 
DE21 7AR (GB)i KING, Steve, P. [GB/GB]; 30 Wil- 
low Drive, Newhall, Swadlincote, Derbyshire DE11 
ONW (GB). KING, Dennis, M. [GB/GB]; 94 Western 


Road, Mickleover, Derby, Derbyshire DE3 5GQ (GB). 
TARASSENKO, Lionel [GB/GB]; 68 Old Road, Head- 
ington, Oxford, Oxfordshire OX3 7LP (GB). HAYTON, 
Paul [GB/GB]; 1 Langford Cottages, London Road, 
Bicester, Oxfordshire OX6 0JN (GB). UTETE, Simukai 
[ZW/GB]; do St. Hugh's College, Oxford, Oxfordshire 
OX2 6LE (GB). 

(74) Agents: ARMITACE, Ian, A. et al.; Mewburn Ellis, York 
House, 23 Kingsway, London, Greater London WC2B 6HP 
(GB). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, FI, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 
MX, MZ, NO, NZ, PL, PT, RO, RU, SD, SE, SG, SI, SK, 
SL, TJ, TM, TR, TT, TZ, UA, UG, US, UZ, VN, YU, ZA, 
ZW. 

(84) Designated States (regional): ARIPO patent (GH, GM, 
\- KE, LS, MW, MZ, SD, SL, SZ, TZ, UG,%£W), Eurasian 
patent (AM, AZ, BY, KG, KZ, MD, RU, Tj|nM), European 

[Continued'vn next page] 


= (54) Title: MONITORING THE HEALTH OF A POWER PLANT 


Vibration 
Amplitude 


== Speed 


^ Next Speed 



3Z 


Change in Amplitude 
at current phase 

Change in phase at same 
Amplitude 


Phase 

^ (57) Abstract: A method for monitoring the health of a system comprises performing at each of a plurality of times the steps of: 
constructing a condition signature (30) from a plurality of condition indicators including (a) a plurality of vibration measurements 

£j acquired from the system or (b) one or more vibration measurements and one or more performance parameter measurements ac- 
quired from the system; predicting a normal signature (32) from a model defining one or more inter-dependencies between said 

^ condition indicators, the normal signature corresponding to the condition signature for a healthy system; comparing the condition 
signature with the normal signature; and registering an event if the condition signature differs from the normal signature by more 

^ than a predetermined threshold. 


wo 02/03041 ai i mil ifflifii d uiiniifli on r d muni umiinnminij iiiiid iiii dii mi 


patent (AT, BE, CH, CY, DE, DK, ES, H, FR, GB, GR, IE, 
IT, LU, MC, NL, PT, SE, TR), OAPI patent (BF, BJ, CF, 
CG, CI, CM, GA, GN, GW, ML, MR, NE, SN, TD, TG). 

Published: 

— with international search report 


— entirely in electronic form (except for this front page) and 
available upon request from the International Bureau 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 


WO 02/03041 


PCT/GB01/03020 


1 

MONITORING THE HEALTH OF A POWER PLANT 

Field of the Invention 

This invention relates to methods and data processing 
5 systems for monitoring the health of a system. The methods 
and data processing systems of the invention are 
particularly, although not necessarily exclusively, 
suitable for monitoring the health of power plant, 
including for example gas turbine, spark ignition and 
10 compression ignition internal combustion engines. 

Background 

The health of a system can be considered a measure of 
the condition of a system against expected norms- A 
T5 healthy, system is one whose condition .closely matches , : 

fy: , expectations, whereas an unhealthy system is orve whose ' - 
conditioh differs from what would be expected; Indicating- 
. .for example deterioration of, ior a. possible, problem with.-: 
the system/,, The ability to monitor' the health of a system 

20 can therefore allow such deterioration and/or problems to 
be detected and, if necessary, addressed at an early stage. 

For example, US-5684718 describes a non-real time 
system for monitoring the operation of an electric 
generator in which vibration and load data are combined to 

25 produce a single signal which is then compared with stored 
data representative of maximum acceptable combinations of 
the two parameters. Essentially the system is an automated 
"look up table'' which issues warnings when vibrations have 
exceeded acceptable limits. 

30 In order to determine the condition, and consequently 

health, of a system, it is normal to monitor and analyse a 
series of measurable indicators which themselves reflect- 
aspects of the condition of the system. For instance, 
talcing the example of a gas turbine, one might monitor 
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performance parameters such as turbine and compressor 
operating temperatures and pressures and spool speeds. To 
obtain a fuller overall picture of the engine's condition, 
these performance parameters can be supplemented with 
5 further condition indicators including, for example, 

vibration measurements and measurements of the particulate 
material entrained in the circulating oil. 

Particularly with complex mechanical systems such as 
gas turbines, the number of indicators that must be 
10 monitored to obtain a useful overall picture of the system's 
condition can be high. This in turn means that the task of 
analysing the complete series of indicators to determine 
the health of the engine is a complex one, typically 
requiring a skilled expert to analyse the data off-line. 
15 ; Taking .again the example of a .gas turbine,: it is known 

,. X • . f br^example : t o collect perf ornpnce^nd, vibration <^ta V vfrOTi 
■ v the engine over time tp be analysed : off- line by one or more 
. experts.;. Typically the performance d^tai will -b^ compared 
with, simulated data for the same engine and, based on this 
20 comparison, an expert will form a view as to the health of 
the engine. Additionally, a small amount of vibration data 
will be reviewed, giving a superficial view of gross 
changes in engine behaviour. If a problem is detected, the 
vibration data may then be analysed in more detail, often 
25 by another expert, to look for any abnormal indications 
which might be symptomatic of underlying mechanical 
problems which could lead to a loss of health and 
operability . 

30 Summary of the Invention 

It is a general aim of the present invention to 
provide methods and data processing systems that facilitate 
the acquisition and analysis of condition indicators in a 
manner such that the overall health of a system can be more 
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readily assessed. 

Accordingly, in general terms a first aspect of the 
invention provides a method for monitoring the health of a 
system, comprising: 
5 constructing a condition signature from a plurality of 

measured condition indicators acquired from the system; 

comparing the condition signature with a normal 
signature, corresponding to the signature for a healthy 
system; and 

10 registering an event if the condition signature 

differs from the normal signature by more than a 
predetermined threshold. 

The term "signature", as used herein, pertains to the 
values of a plurality of condition indicators merged or 

.1-5. .fused into a: unit or quantity , such as a set, vector or 
scailar. In the example of -a r ,ysetor signature, the , 

i ■■■)•■ indicators may correspond to respective elements of the 
vector. ./; In the example of a scalar signature, the 
. magnitude of the. .scalar may ;be .determined by a mathematical 

20 function which acts upon the indicator values. 

By merging or fusing the condition indicators into a 
single signature in this manner, and providing a normal 
signature with which the fused data can be compared, the 
task of assessing the health of a system is greatly 

25 simplified. In particular, since the detection of an event 
amounts to an indication of a potential problem or an 
unhealthy system (i.e. a system condition that differs from 
what would normally be expected) , the monitoring of the 
health can be largely- automated, removing, or at least 

30 minimising, the requirement for expert input during the 
monitoring process. This in turn means that it becomes 
feasible to continuously monitor the health of a system, 
and to provide useful information about the health of the 
system in real time during operation. 
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Preferably the condition indicators that are combined 
to form the system condition signature include operational 
parameters, which in the case of a mechanical system may be 
speeds, pressures (e.g. gas pressures, oil pressures) and 
temperatures for example. Other. useful parameters may 
include what might be conventionally thought of as control 
or status parameters. For convenience, such parameters 
will be referred to using the single label of "performance 
parameters" in the following text. 

Additionally, to obtain a fuller picture of a 
mechanical system's health, it is particularly preferred 
that the signature includes one or more condition 
indicators related to the vibration of the system. 

Put, more generally, the condition indicators from 
which the system, condition signature,.^ 

derived r<m : two ,pr more disparate source s,apf data. >S :^hisv 
illustrates a; ; parti cul<*r strength; of Jbh^s (approach d2n that H 
a ^reat' variety ;o£ different forms of condition indicator ." 
data can be encompassed in -the system- condition^ signature, 
providing a more comprehensive measure of the system's 
health than has previously been possible without multiple 
analyses. 

Preferably at least three condition indicators are 
used to construct the condition signature. More preferably 
at least 10 and even more preferably at least 20 condition 
indicators are used to construct the condition signature. 

In a preferred embodiment, the system comprises a gas 
turbine engine. 

The normal signature for the system can be derived 
from a predefined model of the system that is being 
monitored. This model can itself be developed off-line and 
then fixed for the duration of the operation of the health 
monitoring method. More preferably, however, the model is 
designed to be refined as the method proceeds in order that 
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it might be better tuned to a specific system. 

Whichever approach is adopted, it is particularly 
preferred that the model is a "learnt model" developed 
using a data-driven, or at least partially data-driven 
5 approach. That is to say the learnt model learns from 

training data comprising series of the condition indicators 
which have been labelled as normal (i.e. healthy) or 
abnormal (i.e. unhealthy) as the case may be. In fact, it 
is often the case that normal data is far more readily 
10 available than abnormal data and therefore the training 

data may only include examples of normal data. This still 
results in an effective model, because subsequent events 
can then be identified as departures from the learnt model 
of normality. 

15:** •■ . The normal signature for the healthy system may, be 
ir< predicted .from a; model defining onefcor more;; inter-, s' 
.A, .-.dependencies -between the condition vindicators s^<Thisj. 
fr -enables the model to specify a. continuous boundary: in R*- . v: 
dimensional space (where.; each dimension relates rto';one >of ;N 

20 condition indicators) corresponding to the limits of 

healthy system operation. This is in contrast to "look up 
table" approaches for setting the limits of healthy system 
operation which do not capture the (often complex) inter- 
relationships and correlations between condition 

25 indicators. 

So, for example, it is often the case that the onset 
of a problem or failure in a particular system manifests 
itself in small changes to a number of condition indicators 
which individually, however, remain in their respective 

30 allowable ranges. The "look up table" approach, which is 
only able to sense gross shifts in individual parameters, 
would fail to recognise that a problem or failure had 
occurred. In contrast, when the condition signature for 
the healthy system is predicted from a model defining one 
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or more inter-dependencies between the condition indicators 
the several small changes in the condition indicators may 
have the cumulative effect of driving the condition 
signature outside normal boundary in N-dimensional space. 

Preferably, the predetermined threshold corresponds to 
a statistically significant departure or variance from 
normality as defined by the normal signature. Thus, in the 
example of a normal signature provided by a learnt model, 
further development of the model (e.g. due to the input of 
more training data) will result in a corresponding 
variation in the predetermined threshold. 

In one embodiment, the invention provides a method for 
monitoring the health of a system, which comprises 
performing at each of a plurality of times the steps of: 

> construct of 
^condition. Vindicators includii^g (a) alplprality of vibration; 
measurements Required vfrom the sy^t^ <^v:'(b)^ ^one^or ;more ; ^ 
vibration: measurements . and one or more performance i r, 
parameter, measurements acquired f hom the system; H > 

predicting a normal signature from a model defining 
one or more inter-dependencies between said condition 
indicators, the normal signature corresponding to the 
condition signature for a healthy system; 

comparing the condition signature with the normal 
signature; and 

registering an event if the condition signature 
differs from the normal signature by more than a 
predetermined threshold. 

The model may comprise a matrix (e.g. a covariance 
matrix) with one or more non-zero off-diagonal terms to 
define the inter-dependencies. The step of comparing the 
condition signature with the normal signature may then 
involve calculating a value for the normalised innovations 
squared (NIS) which is defined below in the description of 
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the Embodiments". 

Alternatively, the model may comprise a neural 
network. If there are N condition indicators, one 
embodiment is a neural network which is trained to predict 
5 the value of the N** 1 from the other W-l indicators. The 
step of comparing the condition signature with the normal 
signature may then involve calculating a prediction error 
which is e.g. the square of the difference between the 
predicted value for N and the actual value. There may be N 
10 of these predictive networks operating in parallel for each 
of the condition indicators. In this case the total 
prediction error can be the sum of the prediction errors of 
each of the networks. In another embodiment, a neural 
network is trained to predict a subset of Nl condition 
15 indicators (such v as -vibration values, e.g. at ! a number of 

key frequencies) from another subset 1 of N2 condition 
■ft*. indicators '(such as the performance parameters) , where Nl + 
• • W2 = W. - ^ • ■■■ * VU' ' ■:■ ?v-?. -' < 

Preferably the times define successive intervals of at 
20 most 1 sec duration (i.e. a 1 Hz repetition frequency). 

More preferably the times define successive intervals of at 
most 0.2 sec duration (a 5 Hz repetition frequency), even 
more preferably at most 0.1 sec (a 10 Hz repetition 
frequency) . By acquiring and processing the condition 
25 indicator data at such rates, it is possible for the method 
to monitor the health of the system in real time. 
Therefore, if an abnormal event is registered at any time, 
immediate and appropriate action can be taken by the system 
operator. This can be particularly advantageous for the 
30 operation of safety critical plant such as aero gas turbine 
engines . 

The data acquisition rate can, however, be 
significantly faster than the processing rate. For example 
the data acquisition rate may be in the range 20 Hz to 80 
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kHz. Only a subset of the acquired data may then be 
processed. 

Thus in another embodiment, the invention provides a 
method for monitoring the health of a system, which 
comprises performing at each of a plurality of times 
defining successive intervals of at most 1 sec duration the 
steps of: 

constructing a condition signature from a plurality of 
condition indicators including (a) a plurality of vibration 
measurements acquired from the system or (b) one or more 
vibration measurements and one or more performance 
parameter measurements acquired from the system; 

predicting a normal signature corresponding to the 
condition signature for a healthy system; 

-comparing the condition signature:. with; the normal :- 
•signature;- -and k wJ&-: ' <> "•>:::.£ ■ >* • V^fc.- -;■ ■ 

irfregisteringyan. event if the -condition signature • M r-i. k, 
differs from the normal; signature by more than: ai^j > 
predetermined threshold ^ r > - :s ';«<>.■ j>V,.- 

Where the condition signature is comprised of data 
from disparate sources, for instance performance and 
vibration data, a problem occurs in that the data may well 
not be synchronised in time. If this asynchronous data is 
combined to form the signature, a distorted picture of the 
system's health may well result. For similar reasons, 
training data used to develop a model of normal system 
behaviour should also be synchronised if distortions in the 
model are to be avoided. 

Thus preferably, the condition indicators are 
synchronously acquired from the system to a synchronisation 
imprecision of at most 1 sec. More preferably the 
synchronisation imprecision is at most 0.1, 0.075, 0.0625 
or 0.02 sec. By "synchronisation imprecision" we mean the 
maximum difference between the acquisition times of each 
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pair of condition indicators forming a particular condition 
signature. Desirably, the measurements are acquired from 
the system at a synchronisation imprecision which is less 
than the duration of the successive time intervals, e.g. if 
5 the time intervals are of 0.2 sec duration, the 

synchronisation imprecision may be at most 0.075 sec. 

The invention also provides a data processing system 
for monitoring the health of a system suitable for 
performing the method outlined above. In general terms the 
10 data processing system comprises: 

data acquisition means for acquiring a plurality of 
measured condition indicators from the system; 

processor means for constructing a system condition 
signature from said plurality of measured condition 
15 .indicators; ,-. . f?f . , c - • 'r: ., 

. ; , , comparator means, for comparing the system condition 

v y . signature-. with a predefined .normal signature, corresponding 
to the signature for a healthy system; and ;: 

means for. registering an event if the comparator..- ; 
20 indicates that the system condition signature differs from 
the normal signature by more than a predetermined 
threshold. 

The data processing system may further comprise a 
display means for displaying (a) one or more of the 

25 condition indicators, (b) the result of the comparison of 
the system condition signature with the normal signature 
and/or (c) an alert signal when the comparator indicates 
that the predetermined threshold has been transgressed 
(i.e. an event has been registered). 

30 In one embodiment, the invention provides a data 

processing system for monitoring the health of a system, 
comprising: 

data acquisition means for acquiring a plurality of 
condition indicators from the system at each of a plurality 
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of times, the condition indicators including (a) a 
plurality of vibration measurements or (b) one or more 
vibration measurements and one or more performance 
parameter measurements; 

processor means for constructing a condition signature 
from said condition indicators and for predicting a normal 
signature corresponding to the condition signature for a 
healthy system, the normal signature being predicted by a 
model which defines one or more inter-dependencies between 
said condition indicators; 

comparator means for comparing the condition signature 
with the normal signature; and 

registration means for registering an event if the 
comparator indicates that the condition signature differs 
f rom the ; normal signature by more thanr.a . predetermined 
threshold. , ; • > v .. > ■ _ ■ . ^ : ; „, , : . 

yl , ^n - another embodiment, *h^ • 
processing system for mohito;ring the- health of ;;a system, i 
comprising:. • v--; >;;>.>> . ^ ■> ; . -i 

data acquisition means for acquiring a plurality of 
condition indicators from the system at each of a plurality 
of times defining successive intervals of at most 1 sec 
duration, the condition indicators including (a) a 
plurality of vibration measurements or (b) one or more 
vibration measurements and one or more performance 
parameter measurements; 

processor means for constructing a condition signature 
from said condition indicators and for predicting a normal 
signature corresponding to the condition signature for a 
healthy system; 

comparator means for comparing the condition signature 
with the normal signature; and 

registration means for registering an event if the 
comparator indicates that the condition signature differs 
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from the normal signature by more than a predetermined 
threshold. 

A further aspect of the invention addresses the 
problem of the synchronous acquisition of the condition 
5 indicators. The invention proposes to associate time 

stamps (based on a common clock) with the acquired data and 
to synchronise the data on the basis of these time stamps. 
Accordingly, in this aspect, the invention provides a 
method of synchronising two or more data streams, each data 
10 stream comprising a series of sequentially acquired data 
elements (and relating e.g. to a respective condition 
indicator of the previous aspect), the method comprising: 

associating a time stamp with each data element of 
each stream, the time stamp identifying the time of 
,15 acquisition of the . data element on the basis of a clock ■ ., ■ 
k ■■/ common to' all data streams; .-■ ■ { 

x selecting a first data element from a first stream and 

inspecting its time stamp; 

conducting- a search of; the data elements of .-the or 

20 each other stream to identify the data element in the or 
each other element having an associated time stamp closest 
to that of the selected element of the first stream; and 
marking said identified data element of the or each 
other stream and said selected element of the first stream 

25 as being synchronised with one another. 

Because, the relative acquisition times of the data 
elements are generally more significant than their absolute 
acquisition times, the common clock may be operate within 
an absolute or relative framework. In an absolute 

30 framework one clock provides the time stamp for each data 
element of each data stream. In a relative framework each 
data stream has its own clock, and one of the clocks is 
selected as the reference clock against which the 
acquisition times of the other data streams are measured. 
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It may be convenient to use a mixture of absolute and 
relative frameworks. For example, if the data streams 
relates to performance parameter and vibration 
measurements, the performance parameter measurements may be 
time stamped from one clock and the vibration measurements 
from another clock. 

The process can be repeated until the data elements in 
the first stream have been exhausted. In any subsequent 
processing of the data that is reliant on using 
synchronised data streams, only those data elements marked 
as being synchronised with one another are used. 

In the case where the acquisition rates of the data 
streams differ from one another, it is preferred that the 
first stream, with which the other streams are 
synchronised, is, chosen., to be^ the stream having the lowest 
Requisition rate,. . r; : :i : • ,, : v. .,/.-.;. 1; ■ ^ . . /. ,.•„.,.,....,..• 

'P. ../ .The :. invention .further. prp data probes sing,. 

system for synchronising two r or more data streams, each ; -V 
data stream comprising a, series.; of sequentially acquired-; 
data elements, comprising: 

means for associating a time stamp with each data 
element of each stream, the time stamp identifying the time 
of acquisition of the data element on the basis of a clock 
common to all data streams ; 

means for selecting a first data element from a first 
stream and inspecting its time stamp; 

means for conducting a search of the data elements of 
the or each other stream to identify the data element in 
the or each other element having an associated time stamp 
closest to that of the selected element of the first 
stream; and 

means for marking said identified data element of the 
or each other stream and said selected element of the first 
stream as being synchronised with one another. 
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Brief Description of the Drawings 

The various aspects of the invention will be further 
described by way of example with reference to the 
5 accompanying drawings, in which: 

Fig. 1 schematically illustrates an exemplary data 
structure that can be adopted for operation of the second 
aspect of the invention; 

Fig. 2 shows a neural network architecture for a 
10 learnt model for operation of the first aspect of the 
invention; 

Fig. 3 shows a graph of the prediction error for the 
learnt model on a set of test data corresponding to a 
period of normal operating conditions for an engine; 
415 Fig. 4 shows a graph of. the prediction ;error for the ; 

?■ learnt model for a further engine operating period in which 
'4 - • the engine experienced a bird strike; '■: 

Fig. 5 shows the learning curve' for. a: simple example 
of a system model for operation of the first aspect -of the 
20 invention; 

Fig. 6 shows a comparison of observations and modelled 
estimates for a shaft speed measurement, illustrating 
evolution of the model of Fig. 5. 

Fig. 7 shows the measured low pressure shaft speed 
25 (N1V) for the period of the test data from a more elaborate 
example of the system model, 

Fig. 8 shows the value for the NIS over the same 
period as Fig. 7, 

Fig. 9 shows the values for the 13 condition 
30 indicators and the NIS over the same period from a further 
example of the system model, and 

Fig.. 10 shows a schematic example of an on-the-engine 
health monitoring system. 
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Description of the Embodiments 

The embodiment described below is an example of a data 
processing system employing both aspects of the invention 
discussed above. More specifically, it is a system for 
5 synchronous acquisition, analysis and display of 

performance parameters and vibration data from a power 
plant (e.g. a gas turbine), for monitoring the health of 
the plant. 

In accordance with a preferred form of the second 
10 aspect, the performance and vibration data streams are 
synchronised in real time and, in accordance with a 
preferred aspect of the first aspect of the invention, 
these data are combined or fused to construct a signature 
for the system that can be compared to a signature derived 
15 ■ from a model representing a healthy power, plant , in order- 
f \- to provide; anomaly/event - detection and .hence fault. , .-. - ■ , w ;., 
ii ; : -diagnosis.;, ; . . - :\ •■ : # * v-^, - , 

•v ; • ' "* ? he following discussion focusses on an -application of 
* the system. to monitoring the health of • a gas -turbine . aero- 
20 engine, but it will be appreciated that the methods can be 
adapted to other power plant, including for example ground- 
based and marine gas turbines, and spark ignition and 
compression ignition internal combustion engines, as well 
as other mechanical, thermodynamic, fluid, electrical or 
25 electronic systems. The system acquires performance 

parameters from the gas turbine digitally via an ethernet 
link at a rate between 20 and 40Hz. Typical performance 
parameters are measurements of pressure, temperature, 
thrust, altitude or Mach number. Vibration data is 
30 acquired from analogue vibration transducers which are 
sampled at user-selectable sampling rates (from 625Hz to 
80kHz) via an analogue-to-digital converter. The amplitude 
spectrum of the vibration data is generated using the Fast 
Fourier Transform every 0.2 sec. 
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The performance and vibration data streams are 
asynchronous and stored in separate files together with the 
corresponding timestamps. During review, as data is loaded 
into memory, synchronisation is performed between the 
5 performance and spectrum data on a line by line basis . 

Markers 10,12 (see Fig. 1) are kept which record the last 
synchronised line in the vibration and performance data 
ring buffers 14,16. When new data is available in memory, 
the timestamp for the next vibration spectrum line is 
10 examined. The synchronisation algorithm starts from the 
last previously synchronised location in the performance 
data and searches forwards or backwards based on the 
timestamps of the performance data (accurate to 0.05 sec) 
until the closest matching timestamp in the performance 
lib data ring buffer 16 is . identified. This location in the 
i; performance data is recorded as. 1 being synchronised with the 
■ ; line in-the vibration ring- buffet 14'. : The algorithm then 
proceeds to the next- line -in the vibration ring buffer 14 
(0.2 sec later) arid so on; until there is -no more data 
20 available to synchronise. 

Clearly, therefore, if the performance parameters are 
acquired at 20 Hz (i.e. at 0.05 sec intervals) the 
synchronisation precision is 0.075 sec (i.e. half the 
acquisition interval added to the accuracy of the 
25 timestamps) and if the performance parameters are acquired 
at 40 Hz (i.e. at 0.025 sec intervals) the synchronisation 
precision is 0.0625 sec. 

Considering the synchronisation algorithm in a little 
more detail, it can be seen from Fig. 1 that the algorithm 
30 maintains a synchronisation table 18 that gives the index 
of the performance data entry that matches each vibration 
data line. The algorithm uses variables to mark the latest 
synchronised data in each buffer. The operation of the 
algorithm can be summarised by the following 'pseudo code': 


35 
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15- 


Initialise the latest synchronised markers 
to the start of the vibration and performance 
data. 

2. Loop while there is more data in both ring 
buffers. 

(a) Starting from the latest synchronised 
data item in each ring buffer, examine the time 
stamp, t, on the next entry in the vibration 
ring buffer. 

(b) Search forward in the performance ring 
buffer until a time stamp greater than t is 
found. Select between this entry in the 
performance ring buffer and the previous entry 
for one which .is closest to t andi record the 


>0 


15 


0 


0nGe synchronised^ the analysis of this performance 
and: vibration- data relies on constructing ••models of normal 
jet engine behaviour and then detecting an event or an 
abnormality with respect to these models . 

Traditional aircraft engine monitoring systems are 
based on two distinct processes: the use of vibration 
signatures to indicate engine state, and a separate 
procedure, gas-path analysis, which is employed for 
determination of state from performance parameters. In the 
approach described now, however, performance-related 
parameters such as pressure and temperature can be fused 
with vibration data (such as tracked order vectors - the 
elements of which correspond to the signal amplitude in a 
narrow range centred on the main vibration frequencies for 
each shaft of the turbine) . The aim is to take advantage 
of disparate sources of data to form a more comprehensive 
picture of engine state during normal operation. This in 
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turn should allow a wider range of deviations to be 
identified. 

Furthermore, it is proposed to employ learnt data- 
driven models to model a normal engine. Thus, although 
5 models of the engine system are used, these are not fixed. 
Instead, they evolve with acquired training data. This 
offers the important advantage of robustness. 

The methods of data analysis described below may be 
termed "novelty detection". An advantage of the methods is 
10 that the role of the expert need only be retained in 
classifying training data as abnormal (i.e. novel) or 
normal. The use of Kalman filtering systems in novelty 
detection has been described in e.g. M. Gelb, Applied 
Optimal Estimation. MIT Press 1974. 
...,15.. Two alternative data, analysis methods are described: ... 

below. They - are distinguished by the amount.. ofo prior 
\fc knowledge required to set up the system; In both cases, 1 -r? 
the role of It he expert: heed, onlyr be? retained in ^classifying 
. - .v.;. training, data as novel > or normal. ^ . : v.r- 'v.- ?r - 1 W-> ' 
20 The first method relies on a prior learnt model of 

normality. For example, normal vibration tracked order 
shapes are learnt using a simple clustering model for the 
normal data. The novelty of e.g. the vibration signature 
for an engine under test is assessed by comparing the 
25 closeness of its tracked order signature with the 

prototypical patterns in the clustering model of normality. 
This can be done, for example, by computing the shortest 
. normalised Euclidean distance between the vector encoding 
the tracked order shaped to any of the (prototypical 
30 patterns) cluster centres in the model of normality (see 
Nairac et al, "A System for the Analysis of Jet Engine 
Vibration Data", Integrated Computer-Aided Engineering, 
6(1): 53-65, 1999). If this distance is beyond a previously 
set threshold, the vibration signature as represented by 
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that tracked order is deemed to be outside the bounds of 
normality. In addition to the vibration tracked orders, 
the model of normality for the vibration spectra includes 
the following: sidebands, multiple harmonics, fractional 
harmonics and broadband power. 

The model is illustrated by an example in which a 
neural network having the architecture shown in Fig. 2 was 
developed as the learnt model. 

The neural network had an input layer 30 with four 
nodes for a condition signature consisting of four 
condition indicators measured relating to one shaft of a 
multi-shaft test engine. The condition indicators were the 
vibration amplitude, the phase and the shaft speed all at a 
specified time, and the shaft speed a time increment after 
the specif ied ■ time , . .. . . ,., ^ ... .. . ... . hi , . - 

* The output ;.layer4: 32- of v the, vrig twprk ha*; two nodes ;for 
predicting respectively thei ciiahge^iri, vibr at ion 1 amplitude 
and change in phaseV after the/ time increment -;v - 

The network, had one : hidden layer 34:, each- node ofv 
which contained a Gaussian radial basis function. 

The training phase for network used training data 
obtained from the test engine over a range of normal 
operating conditions. The centres and the spreads of the 
Gaussians were fixed using the cluster analysis described 
above and the weights of the connections between the nodes 
were then iteratively adjusted until the model converged. 

Fig. 3 shows a graph of the prediction error (i.e. the 
sum of the prediction errors of the change in vibration 
amplitude and change in phase) for the model on a set of 
test data which also corresponded to a period of normal 
operating conditions for the engine. This graph provides a 
baseline of prediction error variation against which novel 
events can be judged. 

Fig. 4 shows a graph of the prediction error for a 
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further engine operating period. In this case, however, 
the engine experienced a bird strike. The largest peak in 
the graph corresponds to the moment of bird impact. 
Clearly the model was able to recognise this event. 
5 However, the changed prediction error signal (compared to 
the baseline of Fig. 3) after the event showed that the 
model was also able to detect post-impact abnormal engine 
behaviour. This provides confidence that the model can not 
only detect major events such as bird strikes, but also 

10 more subtle deviations from normality. 

The second method employs a process model which has a 
state vector associated with it (see below) . The 
observation vector (i.e. the condition signature) has 
elements corresponding to measured values of performance 

15 parameters, and. vibration information so, that. two .types of 
data -are.. fused- wdthin the model. . : ; The, fusion of the data is 
performed in real-time with ; a. new output being generated by 
the system ^several times a, second. . . 

_ An important aspect of ..the .use, of this model in , the 

20 system is the use of learning. In a first, off-line, phase 
of learning, a generic model of the engine is learnt. The 
learning is data-driven using an algorithm such as 
Expectation-Maximisation in order to maximise the 
likelihood of the learnt model given the training data. 

25 Once such a generic model has been learnt off-line for a 

particular type of engine, learning can then be applied on- 
line in order to tune the model to an individual engine 
immediately after its pass-off test and after each 
maintenance procedure. Engine deterioration can also be 

30 learnt on-line. The learnt model can be tuned £o different 
flight conditions, such as cruising or landing, in order to 
detect novelty with even more sensitivity and specificity. 

The data-driven learnt model may be integrated with 
existing performance models which rely on the laws of 
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10 


thermodynamics and computational fluid dynamics (knowledge- 
based models) . Such models can therefore be described as 
hybrid models because they are based on the integration of 
learnt and knowledge-based models. 

Looking in more detail at the learnt modelling 
approach, it is based on the application of Expectation 
Maximisation (EM) to parameter estimation in linear 
dynamical systems (see Ghahramani and Hinton, Parameter 
Estimation for Linear Dynamical Systems, Technical Report 
CRG-TR-96-2, University of Toronto, 1996) and to non-linear 
systems (see Roweis and Ghahramani, "A Unifying Review of 
Linear Gaussian Models", Neural Computation, 11, 305-345, 
1999; and Ghahramani and Roweis, "Learning in Nonlinear 
Dynamical Systems Using an EM Algorithm" in Kearns et al. 
15 - (editors) , Advances in Neural Information Processing ^ s- 

''- Systems/ Volume 11, : MIT Press; 1999) ^;/^-*.*^*yj* ••:> 
•"" The ?EM learning algorithm is Applied* to a Kalmah ; - ' •> 

- r- - filter model. - in the linear case;" this is a system with a 
; "measurement process' of the form " ! 3;, vV: «;-^- ; v - 

20 Yd) = Cx(i)+v(i) (1) 

where y(i) is a set of observations of hidden state x(i), c 
is a covariance matrix, and measurement noise v(i) is zero- 
mean and normally distributed with covariance matrix R. 
yd) and x(i) can be the same dimension. Non-zero off- 
25 diagonal terms in C allow the model to account for inter- 
dependencies between the performance parameter and 
vibration measurements of the condition and normal 
signatures. The state equation is 

x(i+l) = Ax(i)+w(i) (2) 
30 with w(i) zero-mean and normally distributed with 
covariance matrix Q. 

At the beginning of the training phase A and C are 
initialised to small random values (e.g. with elements of the 
matrices * 10" s ), and R and Q are initialised e.g. to I. Then 
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during the training phase, for each condition signature y(i) 
in the training set, the method of Roweis and Ghahramani 
applied to equation (2) to derive the most likely values for 
the elements of x(i), and the elements of C, R and v(i) are 
5 iteratively adjusted so that Cx(i)+v(i) converges to the 
respective condition signature (R and Q can be constrained 
throughout to be diagonal matrices) . Convergence can be 
determined by the log-likelihood of the set of observations 
given the model. 

10 Instead of initialising the elements of A to small 

random values, it is also possible to adopt initial values 
that embody existing performance models of engine 
behaviour. After the training phase, the model would then 
be a hybrid of a knowledge-based and a data-driven model. . 
*15 : By fusing these two methods of: -data-analysis, the accuracy 

y • * of prior expert knowledge can be^ combined with 1 ; the • 

| ^ : robustness of data-driven approaches . v : y . ; A -/: v 
When the training phase^hks etiddd and the model is 
; receiving real-time data consisting of 4 a ^sequence' df condition 

20 signatures, the Kalman filter is again used to derive the most 
likely values for the elements of x(i) for each condition 
signature y(i) • However, the elements of C and v(i) are now 
fixed, so Cx (i) +v(i) provides the normal signature for 
comparison with the condition signature. 

25 For example, comparison of the normal signature with 

the condition signature can be on the basis of the 
normalised innovations squared (NIS) . The innovations 
sequence v is the difference between the condition 
signature and the normal signature, so 

30 v(k) = y(k)-c£(k\k-l) (3) 

The innovations should be zero-mean and white. 

The NIS combines the individual innovations sequences, 
NIS(k) - v T (k)S(k)' 1 v(k) (4) 
The individual sequences are weighted by the term S(k)" 1 , 
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the inverse of the innovation covariance given by 
S(k) - C(k)P(k\k-l)C(k) T + E(k) (5) 
where Pftr|]r-i; is the prediction covariance. 

The model: is first illustrated by a. simple example 
(which, does not use vibration measurements) where 
observations are made of the speeds of the three shafts of 
a test engine during cruise. The observed data y is simply 
the state x corrupted by noise, so 

y(i) = xd)-hv(i) (6) 
The observations are used during the learning process, to 
generate a dynamical system model in which A, C, Q and R 
are learned from the data. At the beginning of the 
training phase A and C were initialised to small random 
values and R and a were initialised to J. 

v ' Fig. 5 illustrates the- learning - (log likelihood) plot 
for the system, frig- }6^hows% thfe ^olut i^tf ; of ^estimates of 
shaft T speeds during the: leatning?^ process us in$ the EM • ^ 
algorithm. In the example : shbwri> the ; learning stage lasts 
for the first 25 .Iterations^ H From iteration 25 onwards V 
the system's dynamical properties are determined by the 
learned matrices (which are then kept fixed) . 

Once trained, the systems can be used to detect events 
or abnormalities, that is to say divergences from the 
learnt model of normality. The events of particular 
interest are those that are unexpected, possibly indicating 
a problem with the engine for example. However, 
particularly where the models have been learned only for 
"steady state" parts of the flight envelope (e.g. 
acceleration, cruise and deceleration), transients during 
operation of the engine will also be flagged up as events, 
although they are expected. For example, where a bleed 
valve is opened or closed, the operating condition of the 
engine will exhibit significant differences from a learnt 
model of steady state normality which does not include this 
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event . 

Thus when using such a steady state model, measures 
can be employed to avoid these transient events. For 
instance, since the opening of a bleed valve is an event 
5 that occurs at. a defined point in time, the data collected 
from the engine at that time and slightly either side of it 
(e.g. for 2 seconds before and after) can be eliminated 
from the data analysed by the health monitoring system. 

The approach is next illustrated with a more elaborate 
10 example which uses both performance parameter and vibration 
data. In this example the model was applied to data 
acquired from a test bed-mounted, multi-shaft, aero gas 
turbine engine which was undergoing a cyclic acceleration- 
deceleration test when it experienced an intermediate 
^15 pressure turbine lock-plate,, event . Such . tests are used to 
^ investigate engine behaviour under extreme conditions. The 
f data was acquired asynchronously as described above. ■:■■*< 

A lock-plate : >. event occurred in the engine on test, day 
152-00.: The training-* data was 152-00 data for a period 
20 before the event, and the test data was 152-00 data for the 
period including the event. In this example a 14 
dimensional model (i.e. y(i) and x(i) each had 14 elements) 
was used in which the condition indicator inputs were: 

- The tracked orders of the low pressure (tol) , 

25 intermediate pressure (toi) and high pressure (toh) shafts, 

- The shaft speeds of the low pressure (N1V) , 
intermediate pressure (N2V) and high pressure (N3V) shafts, 

- The ambient (P0V), total inlet (P20V) , high pressure 
compressor delivery (P30V) and exhaust (PEXV) pressures, 

30 - The total inlet (T20V) and high pressure compressor 

delivery (T30V) temperatures, 

- The turbine gas temperature trimmed (TGTTRM) 

- The demanded fuel flow (WFDEM) 

At the beginning of the training phase A and C were 
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initialised to small random values and R and Q were 
initialised to I. 

Fig. 7 shows the measured low pressure shaft speed 
(N1V) for part of the period of the test data, and Fig. 8 
shows the value for the NIS calculated by the trained model 
over the same period. 

The first two sharp troughs in the N1V trace were 
caused by planned consecutive cyclic decelerations. 
Associated with each of these troughs are two NIS peaks. 
These peaks indicated that the engine was not behaving 
normally during the cyclic testing, in fact subsequent 
examination revealed that a lock plate had released earlier 
during the test and as a result abnormal blade rubbing was 
occurring during each of the deceleration cycles. 
15*' During the' third planned deceleration cycle in the ■ ? ; 

;C :^ ibCl covered h v" Fi g s '- ; ? and 8 fife.:/ at approximately data 
Ipbint 12 850) the engine sustained cblade damage which H '•< 
resulted in the sharp NIS peak- arid drop V off ih ' NlV. '--it 
■ ..1: JA.-- ' However, the earlier (all be*it smaller) ; NIS peaks "I" 
20 demonstrate that the monitoring system was able to detect 
the effect of the lock plate release in real time and 
before substantial blade damage was sustained. If such a 
release had occurred in an in-service aero engine, it would 
therefore have been possible to generate an immediate 
25 warning so that timely action (such as engine inspection or 
maintenance) could have been performed. In contrast, 
sudden variations in N1V can occur normally, so N1V alone 
is not a reliable indicator of abnormal behaviour. 

A further example also uses both performance parameter 
and vibration data. Again the model was applied to 
synchronously data acquired from a test bed-mounted, multi- 
shaft, aero gas turbine engine. However, in this case an 
oil seal leak developed in the engine. 

The fault occurred in the engine around data point 
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50410. The training data was from a period before the 
fault, and the test data was for a period including the 
fault. In this example a 13 dimensional model was used in 
which the condition indicator inputs were tol, toi, toh, 
5 N1V, N2V, N3V, P0V, P20V, P30V, PEXV, T20V, TGTTRM, and 
WFDEM. 

Fig. 9 shows the values for the 13 condition 
indicators and the NIS (in the bottom graph) over the 
period including the event. The sharp NIS peak at data 
10 point 50410 again demonstrates that the monitoring system 
was able to detect the moment of the event. Although some 
of the other condition indicators also had peaks at this 
time, by themselves they cannot be reliably associated with 
abnormal (novel) behaviour. 
,15 ,. j - ... >, ..Subsequent ..examination showed that the event was an - 
A : oil; seal leak which allowed engine -oil to enter the 
■:k combustor., Because of .; this leak the engine control system 

reducedi the amount of fuel entering the. combustor (hence 
■ thee drop; in WFDEM ) l to maintain the engine thrust (given by 
20 PEXV/P20) at a constant level. 

After the event, with oil continuing to leak into the 
combustor, the engine was operating outside the scope of 
normal behaviour. This was successfully picked up by the 
NIS trace which after data point 50410 did not return to 
25 its pre-event level. 

Next we consider how a health monitoring system, 
incorporating the model described above, could be installed 
for in-flight analysis of aero gas turbines. 

An on-the-engine system, shown schematically in Fig. 
30 10, could generate of the order of 1 Gb of vibration and 
performance data (consisting mainly of pressures, 
temperatures and shaft speeds) per flight. The vibration 
data is usually analysed in the frequency domain. The 
vibration and performance data, as they are being generated 
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by data acquisition means 20, are temporarily stored in 
ring buffer 22. The data is synchronised and subjected to 
novelty detection in processor and comparator means 24 
which receives a synchronisation signal from data 
acquisition means 20 and the data from ring buffer 22. 
Those sections of the data corresponding to novel events 
are then tagged and recorded with no loss of information 
(i.e. highbandwidth data is recorded) in registation means 
26 which has semi -permanent on-line and/or hard disk 
storage. When the flight is completed the stored data may 
be downloaded and subjected to more intensive ground-based 
analysis. 

The system may also include a display which is driven 
to allow information to be displayed either during 

• acquisition or for review once ; an . acquisition cycle has ...... 

• been completed. ■.. It preferably, includes, the following , . , . 
features: .>>.,, ...»•/_....:,. ^..i ; .., s .. v , ; -. 

- Ability to display the, result of ; the comparison of 
the condition; signature with the^normal - signature, .e.g. in 
the form of the NIS or the prediction error. An unhealthy 
event may be highlighted e.g. with an alert signal. 

- Ability to display a combination of any two of 
vibration spectra, tracked orders, broadband power, 
performance parameters synchronised in time. 

- Ability to extract and plot vibration spectra 
against engine speed. 

- Ability to interrogate and print any of vibration 
spectra, tracked orders, broadband power and performance 
parameters. 

- Automatic detection and display of features from 
vibration spectra (sidebands, harmonics, etc.) 
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CLAIMS 

1. A method for monitoring the health of a system, which 
comprises performing at each of a plurality of times the 

5 steps of: 

constructing a condition signature from a plurality of 
condition indicators including (a) a plurality of vibration 
measurements acquired from the system or (b) one or more 
vibration measurements and one or more performance 
10 parameter measurements acquired from the system; 

predicting a normal signature from a model defining 
one or more inter-dependencies between said condition 
indicators, the normal signature corresponding to the 
condition signature for a healthy system; 
' '15 comparing the 'condition signature with the normal 

v - ,; sigkafcurer arid '•-*- ■ - 7 VT " ' "•• * ' 

'* ; '^registering an event if the condition; signature 
differs f roirr the normal^ signature by more than a 
predetermined threshold. 

20 

2. A method according to claim 1, wherein the model is a 
learnt model. 

3. A method according to claim 1 or 2, wherein the model 
25 comprises a matrix with one or more non-zero off-diagonal 

terms to define said interdependencies. 

4. A method according to claim 3, wherein the step of 
comparing the condition signature with the normal signature 

30 involves calculating a value for the normalised innovations 
squared. 


5. A method according to claim 1 or 2, wherein the model 
comprises a neural network. 
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6. A method according to claim 5, wherein the step of 
comparing the condition signature with the normal signature 
involves calculating a prediction error. 

7. A method according to any one of the previous claims, 
wherein said times define successive intervals of at most 1 
sec duration. 

8. A method for monitoring the health of a system, which 
comprises performing at each of a plurality of times 
defining successive intervals of at most 1 sec duration the 
steps of: 

constructing a condition signature from a plurality of 
condition indicators;, including (a), ,a plurality . of vibration 
measurements^ acquired from $he: system^or (b);-.one .or more \ 
vibration measurement^ or : more • performance ^ •? 

parameter measurements; acquired from the systems ; u 

v - predicting a^nomal signature ^corresponding, to the ■,■ 
condition signature for a healthy system; 

comparing the condition signature with the normal 
signature; and 

registering an event if the condition signature 
differs from the normal signature by more than a 
predetermined threshold. 

9. A method according to claim 8, wherein the normal 
signature is predicted from a model defining one or more 
inter-dependencies between said condition indicators 

10. A method according to claim 9, wherein the model is a 
learnt model. 

11. A method according to claim 9 or 10, wherein the model 
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comprises a matrix with one or more non-zero off-diagonal 
terms to define said interdependencies . 

12. A method according to claim 11, wherein the step of 

5 comparing the condition signature with the normal signature 
involves calculating a value for the normalised innovations 
squared. 

13. A method according to claim 9 or 10 f wherein the model 
10 comprises a neural network. 

14. A method according to claim 13, wherein the step of 
comparing the condition signature with the normal signature 
involves calculating a prediction error. 

. 15. A- method ^according to any one of the previous claims, 
■y.&:.. wherein ,said;, measurements are. synchronously acquired from 
thefi^system^oW synchronisation imprecision of at most 1 



20 

16. , A method according to any one of the previous claims, 
wherein the system comprises a gas turbine engine. 

17. A data processing system for monitoring the health of 
25 a system, comprising: 

data acquisition means for acquiring a plurality of 
condition indicators from the system at each of a plurality 
of times, the condition indicators including (a) a 
plurality of vibration measurements or (b) one or more 
30 vibration measurements and one or more performance 
parameter measurements; 

processor means for constructing a condition signature 
from said condition indicators and for predicting a normal 
signature corresponding to the condition signature for a 
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healthy system, the normal signature being predicted by a 
model which defines one or more inter-dependencies between 
said condition indicators; 

comparator means for. comparing the condition signature 
with the normal signature; and 

registration means for registering an event if the 
comparator indicates that the condition signature differs 
from the normal signature by more than a predetermined 
threshold. 

18. A data processing system for monitoring the health of 
a system, comprising: 

data acquisition means for acquiring a plurality of 
condition indicators from the system at each of a plurality 
of times defining successive intervals of at most 1 sec 
duration,: the ^ 

plurality of - vibration ^measurement or v m0 re 

vibration measurements &hd ' one : or more^ performance 
parameter measurement's ;% -i^y"-'?? ■ ^ ^ - ; - 1 

processor means for constructing a condition signature 
from said condition indicators and for predicting a normal 
signature corresponding to the condition signature for a 
healthy system; 

comparator means for comparing the condition signature 
with the normal signature; and 

registration means for registering an event if the 
comparator indicates that the condition signature differs 
from the normal signature by more than a predetermined 
threshold. 

19. A method of synchronising two or more data streams, 
each data stream comprising a series of sequentially 
acquired data elements, the method comprising: 

associating a time stamp with each data element of 
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each stream, the time stamp identifying the time of 
acquisition of the data element on the basis of a clock 
common to all data streams; 

selecting a first data element from a first stream and 
inspecting its time stamp; 

conducting a search of the data elements of the or 
each other stream to identify the data element in the or 
each other element having an associated time stamp closest 
to that of the selected element of the first stream; and 

marking said identified data element of the or each 
other stream and said selected element of the first stream 
as being synchronised with one another. 

20. A method according to claim 19, wherein the data 
element acquisition rate of the first .stream is lower than 
the acquisition rate of the or, each other streaif. < v.. ; v 

21. A data processing system, for synchronising ; two or more 
data streams, each data stream comprising a series of 
sequentially acquired data elements, comprising ; 

means for associating a time stamp with each data 
element of each stream, the time stamp identifying the time 
of acquisition of the data element on the basis of a clock 
common to all data streams; 

means for selecting a first data element from a first 
stream and inspecting its time stamp; 

means for conducting a search of the data elements of 
the or each other stream to identify the data element in 
the or each other element having an associated time stamp 
closest to that of the selected element of the first 
stream; and 

. means for marking said identified data element of the 
or each other stream and said selected element of the first 
stream as being synchronised with one another. 
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